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EXPERIMENTAL  AND  THEORETICAL  INVESTIGATION 
OF  UNSTEADY  SUPERCAVITATING  HYDROFOILS  OF  FINITE  SPAN 

by 
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Patrick  Leehey 

ABSTRACT 


Measurements  of  the  forces,  moment,  and  cavity  pressure 
on  aspect  ratio  3 and  5 supercavitating  hydrofoils  of 
elliptical  planform  oscillating  in  pitch  have  been  obtained 
at  frequencies  from  5 to  20  Hz.  A region  of  resonant 
behavior  was  found,  which  corresponds  to  a natural  frequency 
of  the  cavity.  A thermodynamic  analysis  of  the  cavity  was 
performed  under  the  assumptions  of  a spherical  cavity  and 
uniform  cavity  interior  properties.  The  calculated  damped 
natural  frequencies  were  compared  with  the  results  taken 
from  peak  oscillatory  lift  experiments.  The  damping  ratio 
of  the  cavity  depends  upon  the  partial  pressure  of  air, 
temperature  and  cavity  size.  A numerical  method  was 
developed  for  a supercavitating  hydrofoil.  Discrete 
vortices  and  sources  are  used  to  formulate  the  equations. 

The  cavity  length  was  iterated  to  get  the  desired  cavitation 
number  over  the  cavitated  planform.  The  numerical  results 
of  an  unsteady  two-dimensional  finite  cavity  were  compared 
with  zero  cavitation  number  theory.  The  problem  of 
supercavitating  hydrofoils  of  finite  span  in  steady  flow 
was  solved  numerically  and  compared  with  analytical  solutions 
and  with  experiments.  Results  of  these  calculations 
indicate  that  this  method  gave  a more  accurate  prediction  of 
lift  cuid  moment  coefficients  than  existing  asymptotic 
theories.  The  numerical  solutions  of  unsteady  hydrofoils 
of  finite  speui  were  also  compared  with  experiments.  The 
amplitude  of  the  forces  was  in  good  agreement,  but  the 
phase  differences  between  theory  and  experiment  was  large. 

The  effect  of  cavity  pressure  variation  should  be  investi- 
gated further. 
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CHAPTER  I INTRODUCTION 

j Considerable  theoretical  and  experimental  work  has 

been  done  on  three-dimensional  supercavitating  hydrofoils, 

1 

jj  but  the  agreement  between  the  theory  and  experiments  is  not 

1 

i fully  satisfactory. 

f 

i Unsteady  two-dimensional  linearized  theories  for 

i 

• zero  cavitation  number  have  been  developed  by  Parkin  (1957) 

and  Woods  (1957) . Leone  (1968)  developed  a two  dimensional 
unsteady  theory  based  on  the  method  of  matched  asymptotic 
expansions,  and  measured  oscillating  lift  for  two-dimensional 
supercavitating  hydrofoils.  A recent  review  of  supercavitat- 
ing flows  around  lifting  hydrofoils  if  that  of  Hsu  (1975). 
Leehey  and  Stellinger  (1975)  obtained  theoretical  predictions 
of  steady  state  forces  and  moments  for  hydrofoils  of  finite 
span  which  compared  well  with  experimental  data.  Widnall 
(1966)  , Tsen  and  Guilband  (1973)  and  Unruh  euid  Bass  (1974) 
studied  linearized  three-dimensional  cavitating  hydrofoils 
using  numerical  lifting  surface  theory.  The  steady  cavity 
length  in  these  studies  is  assumed  to  be  relatively  long  so 
that  the  variation  of  cavity  length  due  to  unsteady  motion 
can  be  neglected  in  the  analysis.  Unsteady  force  character- 
istics were  measured  by  Wetzel  and  Foerster  (1967)  for 
naturally  ventilated  foils  undergoing  either  a sinusoidal 
heaving  motion  or  a harmonic  oscillation  of  a trailing  edge 
flap.  Klose  and  Acosta  (1969)  carried  out  unsteady  force 
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measurements  and  flow  observation  on  superventilated  hydro- 
foils in  heaving  motion.  Their  unsteady  lift  measurements 
showed  appreciable  variation  with  chemges  in  cavity  length. 

The  present  experimental  work  is  concerned  with  non- 
steady motion  of  cavitating  hydrofoils  under  pitch  oscilla- 
tion. A region  of  resonant  behavior  near  the  reduced 
frequency  = 2 was  found,  where  8.  is  the  meem  cavity 

length  measured  from  the  leading  edge  of  the  foil.  For  a 
fixed  oscillation  frequency,  the  eunplitude  of  the  oscillating 
lift  coefficient  peaks  sharply  at  a critical  cavity  length 
(or  volume) . Microphone  measurements  of  cavity  pressure 
show  that  it  is  not  constant  during  the  pitch  oscillation. 

The  resonance  of  the  air-vapor  cavity  has  been  studied 
asssuming  incompressible  flow,  uniform  cavity  interior  proper- 
ties and  a spherical  cavity.  The  liquid  phase  inside  the 
cavity,  viscosity  and  gas  diffusion  have  been  neglected  in 
this  analysis.  A damped  natural  frequency  was  found.  The 
oscillating  eunplitude  depends  on  the  damping  ratio. 

For  long  cavities,  the  oscillatory  forces  are  away 
from  the  region  of  resonant  behavior.  The  unsteady  forces 
are  compared  with  numerical  results,  obtained  from  a lifting 
surface  theory.  The  discrete  vortexes  and  source  method 
is  adopted  to  formulate  the  equations.  An  iteration  scheme 
is  used  to  alter  the  cavity  length  until  the  desired  cavitation 
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nvunber  is  reached  over  the  cavitated  platform.  The  numerical 
cunplitude  of  the  oscillatory  lift  compares  well  with  the 
experiment.  However,  the  phase  difference  between  theory 
and  experiment  cannot  be  neglected. 

Numerical  results  for  lift  and  moment  coefficients  on 
steady  supercavitating  hydrofoils  are  in  good  agreement  with 
steady  flow  experiments  of  Leehey  & Stellinger  (1975)  and 
Maixner  (1977) . The  results  show  that  the  improvement  has 
been  made  after  considering  the  lifting  surface  effects. 
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CHAPTER  II  EXPERIMENTAL  STUDY 
2 . 1 Experimental  Apparatus 

A)  Variable  Pressure  Water  Tunnel 

The  experimental  work  was  conducted  in  the  Recircul- 
ating Variable  Pressure  Water  Tunnel  of  the  Department  of 
Ocean  Engineering,  Massachusetts  Institute  of  Technology. 

Test  section  velocity  is  variable  between  0 and  30  feet 
per  second,  static  pressure  is  variable  between  atmospheric 
pressure  and  3 inch  Hg.  The  test  section  is  20  inches  by 
20  inches  cind  4.5  feet  long.  In  order  to  get  high  oscillat- 
ing forces  and  small  cavitation  numbers,  the  tunnel  velocity 
was  run  as  high  as  practical  { — 30  ft/sec) . The  velocity 
is  indicated  on  a manometer  which  gives  the  differential 
static  head  between  two  locations  in  the  contraction  section. 
This  reading  was  calibrated  by  using  a pitot  tube  in  the 
test  section.  Figure  2.1  presents  an  overall  view  of  the 
working  section.  The  dynamometer  is  mounted  in  the  top 
window  of  the  test  section.  For  these  tests,  two  force  and 
one  moment  component  were  measured. 

B)  Models  and  Shafts 

The  foils  used  in  this  test  were  flat  plate  sections 
with  elliptic  plan  forms,  the  same  used  by  Leehey  and 
Stellinger  (1975),  with  dimensions  as  shown  in  Figure  2.2. 
Each  foil  and  shaft  were  constructed  of  Type  304  stainless 
steel  and  joined  together  by  two  pins. 
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C)  Dynamometer 

The  modified  dynaunometer  was  designed  by  Riva  (1974) 
to  permit  measurements  of  forces  and  moments  on  hydrofoils 
oscillating  in  pitch.  The  pictures  of  this  system  are 
shown  in  Figures  2.1  and  2.3. 

The  shaft  is  supported  horizontally  and  vertically 
by  two  Fafnir  self-aligning  bearings.  These  bearings  are 
mounted  in  a cylindrical-spherical  housing  attached  to 
the  floating  part  of  the  dynamometer,  which  can  be  adjusted 
to  correct  misalignment  of  the  foil  relative  to  the  top 
window. 

The  clamping  and  stopping  system  are  constructed  on 
the  cylindrical-spherical  bearing  housing.  The  clcunping 
system  permits  the  use  of  the  dynamometer  in  steady  experi- 
ments while  the  second  one  prevents  the  foil  from  going  to 
too  large  an  attack  angle  in  case  of  failure  during  unsteady 
tests. 

The  torque  transducer  was  designed  as  an  extension  of 
the  shaft.  It  was  mounted  on  the  top  of  the  shaft  and  with 
set  screws.  The  transducer  is  operated  with  two  foil  strain 
gages  BLH-FAED-06-12-S6 . 

The  oscillatory  movement  of  the  foil  is  produced  by 
two  arms  set  orthoganally  on  2m  imaginary  horizontal  plane, 
and  articulated  by  Fafnir  rod-ends. 

To  provide  lift  and  drag  measurements  without  interfer- 
ence of  one  with  the  other,  the  oscillatory  arms  can  have 
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their  positions  changed,  relative  to  the  center  line  of  the 
tunnel,  which  is  coincident  with  the  center  line  of  the 
dynamometer. 

D)  Drive  Syetem 

The  pitching  apparatus  consists  of  a crank  arm  which 
is  oscillated  sinusoidally  by  means  of  an  eccentric  pin 
arrangement  driven  by  an  electric  motor.  The  motor  speed 
is  controlled  by  a potentiometer. 

A two  level  steel  freune  attached  to  the  flange  and 
side  window  of  the  tunnel  provides  the  necessary  support 
for  the  pulleys,  eccentric  wheel  and  motor.  The  eccentric 
and  two  pulleys  that  drive  the  crank  arm  are  mounted  on  the 
upper  level  of  this  fraune. 

The  pulleys  have  speed  ratio  2 and  3 relative  to  the 
motor.  The  motor  is  a Master  DC  motor  with  1/2  HP  output 
at  2000  RPM.  The  eccentric  is  designed  to  provide  in  the 
following  aunplitudes  of  dynamic  angle  of  attack  of  the  foil: 
1.5°,  2.0°,  2.5°,  3.0°,  3.5°,  and  4.0°. 


2.2  Hydrofoil  Test  Procedure 


I Dyneimic  calibration  of  the  force  gages  were  carried 

out  with  the  aspect  ratio  5 foil.  The  dynamometer,  installed 
in  the  water  tunnel,  was  subjected  to  a sinusoidally  oscilla- 
ting force,  generated  by  the  B & K shaker.  The  force  signals 
from  the  force  gages  and  the  shaker  were  analyzed  with  the 
Spectral  Dynamics  Analyzer,  a heterodyning  filter  with  5 Hz 
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bemdwidth.  The  measurements  covered  the  frequency  range 
from  5 Hz  up  to  20  Hz.  The  force  ratio  and  phase  difference 
between  the  load  cells  emd  shaker  show  that  the  natural 
frequency  of  the  dynamometer  is  above  the  test  range.  It  is 
reasonable  to  represent  the  real  system  by  the  model  with 
mass-spring-dashpot . Damping  is  small  in  this  mechanical 
system  so  that  the  exact  value  is  vinimportant . Details  of 
this  calibration  are  shown  in  Appendix  A. 

A piezoelectric  pressure  transducer  was  mounted  at 
the  tunnel  top  window  about  5"  from  midchord  to  measure  the 
dynamic  pressure  inside  the  cavity.  The  diaphgram  of  the 
pressure  transducer  is  at  1"  below  the  root  of  the  foil. 

An  accelerometer  was  mounted  on  the  crank  arm  at  7"  from  the 
foil  shaft  to  count  the  frequency  and  to  provide  phase 
difference  relative  to  force,  moment  and  cavity  pressure. 
Signals  output  of  transducers  were  amplified  and  then  analyzed 
by  the  Spectral  Dyneimics  Analyzer.  A general  arrangement  of 
the  readout  system  is  shown  in  Figure  2.4. 

The  dyncimometer  was  aligned  at  the  flow  direction  to 
read  lift  and  drag  directly.  Geometrical  angle  of  attack 
was  set  by  aligning  the  flat  side  of  each  foil  with  a wedge 
which  against  the  test  section  side  wall.  Amplitude  of 
oscillation  can  be  set  by  selecting  the  wheel  eccentricity. 
Oscillation  frequency  was  adjustable  from  5 Hz  to  20  Hz  by 
changing  the  resistance  in  the  field  windings  of  the  DC  drive 


motor. 
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Static  pressure  in  the  test  section  was  read  from  a 
mercury  manometer  connected  to  taps  in  the  tunnel  wall  at 
the  beginning  and  end  of  the  test  section,  on  the  center 
line,  and  corrected  for  hydrostatic  head  to  give  static 
pressure  at  hydrofoil  centroid.  Static  pressure  in  the  cavity 
was  read  from  a mercury  manometer  connected  to  a tap  on  the 
suction  side  of  the  foil.  Test  section  static  pressure  was 
variable  between  3"  Hg  and  9.5"  Hg  absolute  pressure. 

Measurements  have  been  made  for  the  forces  and  moments 
on  aspect  ratio  3 and  5 hydrofoils  oscillating  in  pitch. 
Readings  were  taken  of  lift,  drag,  moment,  static  pressure, 
cavity  pressure,  frequency,  oscillatory  forces,  phase  differ- 
ence, water  temperature  and  velocity.  At  each  oscillation 
frequency,  the  static  pressure  of  the  test  section  was 
reduced  in  increments  of  0.6"  Hg  in  the  supercavitating  case. 
Cavity  length  was  photographed  and  confirmed  with  Leehey  & 
Stellinger  (1975) , and  then  the  mean  cavity  length  is  taken 
emd  used  as  the  parameter  in  the  reduced  frequency. 

The  cavitation  number,  , is  defined  as 


- t-fc 


(■2-0 


where  is  the  pressure  for  upstream  at  the  elevation  of 
the  foil  centroid,  is  the  cavity  pressure.  The  cavitation 
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number,  0^  , based  on  vapor  pressure  is  defined  as 


O = — 

''  T?li^  . 


C2-2) 


The  vapor  pressure,  , was  calculated  as  a function  of 

water  temperature  according  to  the  formulae  of  Smith,  Keyes 
and  Gerry  (1934)  . 

In  this  experiment,  the  cavity  pressure  is  taken  from 
the  pressure  tap  on  the  foil. 


2.3  Experimental  Results  & Discussions 

The  comparisons  of  measured  lift  and  moment  for  unsteady 
fully  wetted  flow  with  analytic  solutions  are  shown  in  Figures 
2.5  and  2.6.  The  foils  used  in  the  analytic  solutions  are 
flat  plates  which  differ  fro.,i  the  foils  we  used  in  our  test. 
However,  the  agreement  in  both  phase  and  amplitude  is  good 
for  lift.  The  torsional  stiffness  of  the  dynamometer  was 
changed  when  the  tunnel  pressure  was  reduced.  In  particular, 
the  stiffness  of  the  rubber  seal  which  clamps  the  shaft 
depends  upon  the  pressure  difference  between  the  atmosphere 
and  tunnel  pressure.  Therefore,  the  properties  of  the  rubber 
seal  have  been  measured  for  different  pressure  differences 
across  it,  permitting  the  data  reduction  to  obtain  the  moment 
coefficient  based  on  the  calibration  of  the  torque  transducer 
under  operational  tunnel  pressures. 


Figures  2.7  and  2.8  present  the  oscillating  lift  force 
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and  the  phase  difference  between  the  force  and  displacement. 
There  is  evidence  of  a cavity  resonance  at  reduced  frequency 
co£ /{Ja,^  2.  • Resonance  frequency  decreases  as  the  cavity 
length  increases.  For  longer  cavities,  the  effect  of  this 
kind  of  resonance  is  weak,  because  the  damping  of  cavity 
oscillation  is  greater  for  long  cavities.  An  analysis  of 
this  phenomena  appears  in  the  following  chapter.  From  the 
observation  of  high-speed  motion  pictures  it  appears  that 
the  cavity  length  changed  a lot  at  the  resonance  frequency. 
The  variation  of  cavity  length  depends  on  the  oscillating 
frequency  emd  flow  velocity.  The  phase  difference  between 
lift  force  and  angular  displacement  tends  to  decrease  as 
the  oscillation  approaches  the  resonance. 

Figures  2.9  through  2.12  show  the  moment  coefficients,  drag 
coefficients  and  phase  differences.  The  oscillatory  coeffi- 
cients depend  on  the  cavity  length.  The  peak  variations 
occur  near  the  resonant  frequency. 


Measurements  of  cavity  pressure  show  that  the  cavity 
pressure  is  not  constant  during  the  pitch  oscillation.  The 
amplitude  of  cavity  pressure  oscillation  depends  on  the 
cavity  length  and  oscillation  frequency,  being  greatest  at 
high  frequencies  and  short  cavity  lengths.  The  relationship 
between  unsteady  cavity  pressure  and  oscillatory  lift  cieffi- 
cient  at  two  different  reduced  frequencies  for  each  foil  is 


( 
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CHAPTER  III  ANALYSIS  OF  BUBBLE  RESONANCE 


As  mentioned  in  Chapter  II,  there  is  evidence  of  a 
cavity  resonance  near  the  reduced  frequency  at  2.  Determin- 
ation of  the  natural  frequency  of  the  cavity  flow  system  is 
very  complicated  matters.  For  incompressible  flow,  the 
complexity  of  the  problem  is  greatly  reduced  by  the  assump- 
tions of  spherical  symmetry  and  uniform  cavity  properties. 
The  viscosity  is  neglected  in  the  following  discussion. 
Several  cases  of  the  bubble  dynamics  have  been  discussed  by 
Hsieh  (1965) . Evaporation  is  an  important  process  in  cavity 
oscillation.  The  variation  of  cavity  pressure  in  time  can 
not  be  neglected  as  it  is  in  the  growth  of  small  cavities. 
The  thermodynamic  effects  are  essential  in  the  cavity 
oscillation.  The  simplified  linear  governing  equations  of 
the  cavity  oscillation  about  the  mean  quantities  are 
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The  interfacial  boundary  conditions  (r  = R)  are 


T = T . 
V = R , 

f 


(3-2) 


and 


i if jf.'tc'iJ+r.'jl . 

In  the  foregoing  set  of  equations,  V is  velocity,  p is 
pressure,  T is  temperature,  f is  density,  is  specific 

heat  at  constant  volume,  ^ is  the  coefficient  of  thermal 
conductivity.  Primed  symbols  with  subscript  for  pertur- 
bation quantities,  and  with  subscript  (o)  for  the 
equalibrium  quantities  within  the  cavity,  curled  symbol  for 
perturbation  quantities  outside  the  bubble,  subscript  (o) 
for  equilibrixim  values.  In  the  last  boundary  condition 
equation  we  neglect  the  heat  conduction  through  the  liquid 
phase  within  the  cavity. 

When  the  cavity  is  undergoing  harmonic  oscillation, 
the  solutions  of  Eq.  (3-1)  subject  to  the  boundary  condition 


Eq.  (3-2)  will  be 


V = 


* 


(3-4) 
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where  K is  the  modified  Bessel  function  and  D ( s ) 

is  the  thermal  diffusivity  of  water.  The  velocity  and 
temperature  gradient  at  the  cavity  boundary  are 


dT  _ 


(^-6) 


Since  D is  a small  number  ( c:  i.s  x lo^-f^Vsec.  ) the  argument  of 
the  Bessel  function  is  large.  The  temperature  gradient  is 


approximately 


= -t'(-^)  e 


Vz  iirU 
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Substituting  (3-5) , (3-7)  into  (3-3) , we  get  the  condition 


P.toR, 


(e:L  f.')  r '+ ( if.  K.c; « v ^ c 


+ - 0, 
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If  the  gases  within  the  cavity  satisfy  the  ideal  gas 
equation  fv  = (^T  and  we  neglect  the  diffusion  of  gases 
(i.e.  = constant)  , 

"P.  = P.iR.T. 

After  a simple  calculation,  and  using  (3-5),  one  can  get  the 
variation  of  pressure 
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We  assume  that  the  density  of  the  vapor  corresponds  to  its 
temperature  at  the  given  pressure  according  to  the  equation 


of  state. 


\ = P„iR.T 


The  variation  of  vapor  pressure  is 


t;= 


ii-lo) 


We  assume  the  mass  of  gases  inside  the  cavity  is  constcuit 


PaV  — Gjfist&.nt'. 


From  this  equation  cuid  Eq.  (3.5)  we  get 

P ’ = -to ' 

then 


Substitute  the  last  equation  into  Eq.  (3-10) 


fv  =fv(RvT'+(^T(e'-^or 


The  total  variation  of  cavity  pressure  is  the  summation 
of  the  variation  of  partial  pressures. 

f '=  t;  - ?: , 


where 


r:  ^p.'«T, 


j The  vapor  within  the  cavity  is  not  a superheated  vapor, 

i ; 

f|  A relationship  between  vapor  pressure  and  temperature  was  I 

I I 

I taken  from  the  paper  by  Smith,  Keyes  and  Gerry  (1934) 
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I 1 
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We  can  get  the  equation  of  variation  of  pressure  and  temp- 
erature by  taking  the  derivative  of  pressure  with  respect 
to  temperature 


T,'  = 0, 


where 


(Tc +(c('Tc  - »')  k -x‘]*  J 


Combine  Eqs.  (3-9)  and(3-12) 


r»  R.  to 


(3-/3) 


Equations  (3-8),  (3-11)  and(3-13)  constitute  the  eigenvalue 
problem.  The  natural  frequency  of  the  system  is  determined 
by 


r _l_  3fa,tf(vT  \ 

^ ' f.R.V  ^ 
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or 


Au;*  - d = j r ^ ^ 


O' 14) 


where 


LRtC  (Ry  J^>C>  _ ^ 

“ if/  « 3 3T: 


g ^ ^cfK-  e>;(g.T)  T-.c; 


t,T.  /?.  ■ ^ 

, (RyC  p;l+  v:  l?^ 
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The  typical  natural  frequencies  for  a given  T,  'fao, 

R^,  are  shown  in  Figure  3.1.  We  assume  the  gases  in  the 
cavity  are  air.  A deunped  natural  frequency  was  found. 

-O  * 

At  a fixed  fao,  the  natural  frequency  decreases  when  the 
cavity  volume  increases.  The  air  content  not  only  Increases 
the  natural  frequency  but  also  reduces  the  damping  ratio. 

For  the  damped,  single-degree  of  freedom  case,  the  ratio  of 
deunping  to  critical  deunping  is  equal  to  the  ratio  of  A/i  to  Ahi, 
where  cJi  is  the  imaginary  part  of  natural  frequency,  and  <dn 
is  the  undamped  natural  frequency.  The  ratio  of  oJi  to  cJn  will 
give  some  idea  of  the  damping  for  different  cavities  as  shown 


rr 
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in  Figure  3.2.  This  analysis  will  not  only  tell  us  the 
natural  frequency  but  will  explain  the  qualitativety  oscil- 
lating cunplitude  of  lift  coefficient  shown  in  Figures  2.7 
eind  2.8.  The  cavity  pressure  is  always  higher  than  the 
vapor  pressure.  This  result  is  to  be  expected  as  the  gases 
from  the  water  contribute  to  the  pressure  within  the  cavity. 

The  discrepancy  between  cavity  pressure  and  vapor  pressure 
decreases  with  decreasing  cavitation  number.  Experimental 
results  of  cavity  pressure  measurement  are  shown  in  figure 
4.10.  Further  discussions  appear  in  the  next  chapter. 

Since  the  short  cavaties  have  a higher  partial  pressure  due 
to  non-water-vapor  gases,  damping  ratio  (damping/critical 
damping)  is  small  in  short  cavities.  The  stiffness  of  cavity 
increases  with  increasing  air  pressure  inside  the  cavity,  so 
we  expect  that  they  have  higher  natural  frequencies  and  less 
deunping  ratio.  However,  the  reason  for  shorter  cavities  having  a 
higher  gas  content  is  not  clear. 

The  cavity  volumes  are  calculated  from  the  numerical 
lifting  surface  method.  Under  the  assumption  of  spherical 
voltime  emd  neglecting  the  surface  tension  effect,  the 
natural  frequency  of  the  cavity  calculated  from  Equation  (3-14) , 
and  the  resonance  frequency  under  test  are  shown  in  Figure 
(3.2).  In  the  present  experiment,  only  a few  points  can  be 
plotted  in  the  graph.  The  "loop"  trajectories  in  Klose  & 

Acosta  (1969)  seems  also  associated  with  the  natural  frequency 
of  the  cavity. 
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CHAPTER  IV  NUMERICAL  ANALYSIS 

4.1  Linearized  Theory  for  Three-Dimensional  Supercavitating 
Hydrofoils 

The  linearized  problem  for  a cavitating  hydrofoil  of 
finite  span  in  unsteady  flow  cam  be  represented  by  a set  of 
coupled  integral  equations.  The  derivation  proceeds  either 
from  Green's  theorem  or  from  the  solution  for  the  velocity 
field  induced  by  a distribution  of  vortices  and  sources. 

This  distribution  must  be  chosen  to  satisfy  the  relevant 
boundary  conditions.  The  numerical  model  we  are  about  to 
describe  is  shown  in  Figure  4.1 

The  source  distribution,  q (x,z,t  ),  represents  the 
slope  difference  between  the  upper  cavity  surface  and  the 
mean  camber  line  of  the  foil  or  the  lower  cavity  surface  at 
the  point  (x,o,z).  The  vortex  distributions  (y  (x,z,t), 
6(x,z,t  ))  create  the  perturbation  velocity  differences  in 
the  X and  z components  between  the  upper  and  lower  surfaces 
at  the  point  (x,o,z).  This  leads  to  the  result 

vcx.+o.j.*)  - J,t))  , 

(T  ( u U (4--I) 


The  boundary  conditions  specified  are  the  normal 
velocity  on  the  wetted  surface  of  the  foil,  a constant 
pressure  in  the  cavity,  and  closure  of  the  cavity. 

+0.  i,t)  , (4-2) 

0 = JJ  , 

5t 

where  h(x,z,t)  is  the  instantaneous  y-ordinate  of  the  wetted 
surface.  is  the  surface  of  the  cavity  as  projected  onto 

the  x-z  plcine.  o is  the  cavitation  number.  P and  P are 
free-stream  pressure  and  cavity  pressure,  respectively. 

Geurst  (1961) proved  that  for  two-dimensional  steady  flow, 
the  re-entrant  jet  and  Riabouchinsky  models  for  cavity  term- 
ination reduce  upon  linearization  to  a statement  that  the 
cavity  is  closed  at  its  end.  The  requirement  that  the  source 
strength  be  zero  beyond  the  cavity  trailing  edge  and  the  sum 
of  sources  is  zero  are  necessary  to  meet  this  condition. 

Based  on  the  boundary  conditions,  the  integral 
equations  for  three-dimensional  unsteady  flow  can  be  written 
as  (The  details  of  this  derivation  are  shown  in  Appendix  B) 
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where  S is  the  projection  of  the  foil  surface  on  to  the  x-z 

Sl 

plane  and  is  the  projection  of  the  wake  region  on  to  the 
same  plane.  Note  that  the  projected  cavity  area,  S^,  is 
also  time  dependent.  The  unknowns  are  the  source  strength, 
vortex  strength  and  cavity  surface  location. 

For  steady  flow  at  a given  cavitation  number,  this  set 
of  equations  simplifies  to 


Xr(^-c?-S)  = -7 
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4.2  Linearized  Theory  for  Two-Dimensional  Supercavatating 
Hydrofoils 

In  the  two-dimensional  theoretical  model,  with  an 
incompressible  fluid,  infinite  pressures  are  necessary  at 
f->oo  to  generate  a volume  change  in  the  cavity.  However, 
Benjamin  (1964)  pointed  out  that  any  real  flow  is  necessarily 
three-dimensional  'in  the  large';  thus  the  pressure  singular- 
ity is  an  artifact  of  the  two-dimensional  case.  Based  on 
his  discussion,  the  changing  cavity  volume  of  two-dimensional 
foil  is  introduced  by  the  three-dimensional  outer  solutions. 
The  overall  solution  can  be  obtained  by  matching  the  inner 
and  outer  solutions.  Theoretically,  the  two-dimensional 
unsteady  cavity  can  be  treated  as  the  inner  solution  of  the 
three  dimensional  case.  This  inner  solution  does  not  match 
the  boundary  condition  at  infinity.  The  numerical  method 
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for  unsteady  flow  appears  in  this  section  only  for  the  inner 
flow  field. 

The  development  of  this  problem  proceeds  from  the 
velocity  or  pressure  field  induced  by  a two-dimensional 
distribution  of  vortices  and  sources.  This  velocity  or 
pressure  field  will  satisfy  the  boundary  conditions  which 
require  a given  normal  velocity  on  the  wetted  surface,  con- 
stant cavity  pressure  and  closure  of  the  cavity. 

The  derivation  of  these  integral  equations  is  the  same 
as  for  three-dimensional  cavity  flow,  except  two-dimensional 
sources  and  vortices  are  used.  The  potential  due  to  two- 
dimensional  source  distributed  over  the  cavity  region  is 


t 
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The  two-dimensional  vortex  potential  will  be 


(4-y) 


(4-/0) 


The  integral  equations  which  satisfy  the  boundary 
conditions  are  converted  to 
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These  equations  are  sufficient  to  obtain  the  cavity 
length,  the  source  distribution  and  the  vortex  distribution. 


I 


4.3  Numerical  Method 

The  solution  of  the  problem  is  obtained  by  reducing 
the  coupled  integral  equations  to  a set  of  simultaneous 
algebraic  equations.  Since  the  method  used  in  two-dimensional 
flow  and  three-dimensional  steady  flow  is  a simplification 
of  the  three-dimensional  unsteady  problem,  the  approach  for 
the  latter  case  will  be  discussed  here. 

The  surface  is  divided  into  small  elements.  The  element 
used  in  this  model  is  one  obtained  by  dividing  the  foil  and 
cavity  semispan  into  strips  with  cosine  spacing,  while  the 
chord  is  divided  into  strips  of  constant  spacing  on  the  foil 
emd  on  the  cavity  behind  the  foil.  In  order  to  increase  the 
rate  of  convergence  for  the  nxunber  of  elements  used,  the 
first  two  elements  of  each  chordwise  strip,  near  the  leading 
edge,  have  only  one-half  the  chord  spacing  of  the  rest  of 
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the  elements  on  the  foil. 

The  model  uses  discrete  vortices  and  sources  to 
represent  the  foil  and  cavity.  Each  element  contains  a 
bound  line  vortex,  a trailing  vortex  and  a line  source. 

The  quarter-chord  line  of  each  element  contains  the  bound 
vortex,  and  the  induced  velocities  are  calculated  for  all 
elements  at  their  midspan,  three-quarter  chord  positions. 
James  (1972)  showed  that  for  the  best  efficiency  the  vortex 
and  control  points  should  be  arranged  according  to  the 
Pistolesi  approximation,  i.e.,  at  the  quarter  and  three- 
quarter  chord  points  on  each  element.  The  concentrated  line 
source  is  taken  to  be  a constant  distribution  across  each 
element  at  its  three-quarter  chord  position.  Since  the 
ijource  distribution  is  singular  at  the  cavity  leading  edge 
and  termination  point,  these  control  points  should  be  placed 
away  from  the  singularities.  Therefore,  the  first  control 
point  is  located  at  the  three-quarter  chord  and  the  source 
is  located  at  the  quarter  chord.  The  local  vortex  or  source 
strength  divided  by  the  element  width. 

Since  the  cavity  length  is  unknown,  the  global  closure 
condition  is  not  used  in  the  calculation.  If  it  is,  the 
variation  of  cavity  length  along  the  span  cannot  be  iterated. 
This  difficulty  is  overcome  by  applying  the  closure  condition 


to  individual  spanwise  strips.  This  approximation  is  exact 
for  a symmetric  foil.  The  source  integration  along  a strip 
between  z = - C and  z = C is 
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the  elements  on  the  foil. 

The  model  uses  discrete  vortices  and  sources  to 
represent  the  foil  and  cavity.  Each  element  contains  a 
bound  line  vortex,  a trailing  vortex  and  a line  source. 

The  quarter-chord  line  of  each  element  contains  the  bound 
vortex,  and  the  induced  velocities  are  calculated  for  all 
elements  at  their  midspan,  three-quarter  chord  positions. 
James  (1972)  showed  that  for  the  best  efficiency  the  vortex 
and  control  points  should  be  arranged  according  to  the 
Pistolesi  approximation,  i.e.,  at  the  quarter  and  three- 
quarter  chord  points  on  each  element.  The  concentrated  line 
source  is  taken  to  be  a constant  distribution  across  each 
element  at  its  three-quarter  chord  position.  Since  the 
source  distribution  is  singular  at  the  cavity  leading  edge 
and  termination  point,  these  control  points  should  be  placed 
away  from  the  singularities.  Therefore,  the  first  control 
point  is  located  at  the  three-quarter  chord  and  the  source 
is  located  at  the  quarter  chord.  The  local  vortex  or  source 
strength  is  given  by  the  discrete  vortex  or  source  strength 
divided  by  the  element  width. 

Since  the  cavity  length  is  unknown,  a closure  condition 
is  applied  to  individual  spanwise  strips  in  the  calculation. 
Otherwise,  the  iteration  technique  does  not  converge.  This 
condition  satisfies  the  overall  global  condition.  It  also 
allows  the  variation  of  cavity  length  along  the  span  to  vary. 
This  method  is  exact  for  a symmetric  foil.  The  source  inte- 
gration along  a strip  between  z = and  z = Zj^  is 


Jj 
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where  ( C ) and  5-  (C,t)  are  the  spanwlse  profiles  of  the 
leading  edge  and  cavity  termination  points,  respectively. 
Let  g^  (x,z,t)  represent  the  profile  of  the  cavity-foil 
combination.  Then  the  stripwise  integration  becomes 

f r /■  \ Jt 
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= 0. 
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The  first  integral  is  zero  if  the  cavity  is  closed 
both  at  its  leading  edge  and  trailing  edge.  The  second 
integral  is  zero  due  to  the  symmetry  of  the  cavity.  This 
integral  therefore  represents  the  stripwise  closure  condition. 

The  effect  of  all  singularities  must  be  calculated  at 
each  control  points  (Appendices  C and  D)  . The  Biot-Savart ' s 
law  is  used  to  calculate  the  induced  velocity  at  a control 
point  due  to  a constant  strength  discrete  vortex  segment. 

The  Kirchhoff's  law  and  Kelvin's  theorem  must  be  satisfied 
for  the  vortex  distribution.  Consider  the  vortex  segment 
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(Xj^.o.Zj^)  to  (x2,o,z^)  with  strength  y per  unit  length. 
Defining  (^,t)  as  the  coordinates  of  a general  point  on 
the  vortex,  and  (x,2)  as  the  coordinates  of  a control  point, 
the  induced  velocity  at  (x,z)  due  to  the  bound  vorticity. 


■+■  C0S&  ) ^ 
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The  induced  velocity  due  to  the  trailing  vortex  (x,o,Zj^)  to 
(x3,o,Z2)  with  strength  is 


4-TTir(^0,J,f) 
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The  total  velocity  induced  by  a vortex  distribution 
is  the  sum  of  the  influence  of  all  the  spanwise  and  chord- 
wise  vortices. 
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According  to  the  Bernoulli  equation,  the  pressure 
induced  by  the  source  component  contains  two  terms;  the 
space  derivative  and  the  time  derivative. 


(4-<9) 


The  Pg  due  to  the  line  source  ((Xj^,o,z^)  - can 

be  calculated  by  line  integration 
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where  a,  d.  t.  Y.  and  3 are  defined  the  same  as  in  Eo.  (4-17). 
({•p  and  4*2  are  the  ancles  shown  in  Figure  4.2. 

Eauations  (4-16)  and  (4-20)  are  not  suitable  for 
numerical  computation  if  d//a  becomes  small.  Approximate 
formulas  are  used  for  this  special  case. 

An  iteration  is  used  to  alter  the  cavity  length  until 
the  desired  cavitation  number  is  reached  over  the  cavitated 
planform.  The  steady  flow  solution  is  determined  first, 
then  the  unsteady  problem  is  solved  by  a step-by-step 
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procedure.  Since  the  variation  of  the  cavity  length  is 
unknown,  the  use  of  the  frequency  domain  is  not  suitable. 

At  each  time  step,  a cavity  length  is  assumed  and  the  matrix  . 

coefficients  are  calculated.  Solving  the  matrix  with  approp- 
iate  boundary  conditions,  we  can  get  the  distribution  of 
vortices  and  sources,  and  the  cavitation  number  on  each  strip. 

Then  a new  cavity  length  is  chosen  using  the  calculated  ; 

cavitation  number.  The  length  is  initially  adjusted  by  one 

percent  of  the  steady  state  cavity  length.  Then  either 

extrapolation  or  Muller's  iteration  scheme  (Kristiansen  1963) 

is  used  to  calculate  cavity  length. 

The  linear  algebraic  equations  which  satisfy  the 
boundary  conditions  are  solved  by  using  Gaussin  elimination 
with  equilibration  and  partial  pivoting  (e.g.  Isaacson  and 
Keller  (1966)).  The  convergence  test,  both  in  space  and  time, 
has  been  made.  Tedile  4.1  gives  the  number  of  elements  chosen 

in  the  numerical  calculation. 

4.4  Comparison  Between  Theory  and  Experiment 

Present  numerical  results  for  two-dimensional  super- 
cavitating  flow  are  compared  with  Geurst's  analytic  solution 
(1960)  and  Golden's  numerical  scheme  (1975).  Golden  used  a 
uniform  distribution  of  the  singularity  in  each  element  and 
the  best  choice  of  control  points.  Figure  4.3  simply  confirms 
that  these  methods  of  solution  lead  to  identical  results. 
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The  lift  and  moment  calculation  for  supercavitating 
hydrofoils  of  elliptic  planform  was  performed  and  compared 
with  analytic  solutions  by  Leehey  (1973) , experiments  by 
Leehey  and  Stellinger  (1975)  and  experiments  by  Maixner 
(1977).  Leehey  utilized  the  method  of  matched  asymptotic- 
expansions.  The  theory  is  valid  to  first  order  in  angle 
of  attac)c  and  second  order  in  the  reciprocal  of  aspect  ratio. 
The  experiments  of  Leehey  and  Stellinger  measured  forces, 
moments  and  cavity  length  for  aspect  ratio  3 and  5 super- 
cavitating hydrofoils  of  elliptic  planform.  Maixner 
investigated  the  water  tunnel  wall  effects  on  supercavitat- 
ing hydrofoils  of  finite  span  by  using  a geometrically 
similar  family  of  three  hydrofoils  (AR  = 5) . He  showed 
that  the  previous  tests  by  Leehey  and  Stellinger  were  reli- 
able. The  standard  wind  tunnel  wall  correction  is  adequate 
when  the  foil  span  is  equal  to  or  less  than  half  of  the 
tunnel  depth.  It  should  be  noted  that  the  data  of  Leehey 
and  Stellinger  was  based  on  the  cavitation  number  calculated 
by  using  the  vapor  pressure  rather  than  the  measured  cavity 
pressure.  In  the  comparison  with  present  results,  for  the 
aspect  ratio  5 foil,  experimental  data  for  the  medium  foil 
is  ta)cen  from  Maixner' s results  where  measured  cavity  pres- 
sure were  used. 

Figures  4.4  cuid  4.5  show  the  lift  coefficient,  C^/ot. 
versus  0/o(  for  the  aspect  ratio  5 and  3 hydrofoils,  respect- 
ively. The  agreement  between  theory  and  experiment  is  good  . 
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Por  small  values  of  O/oc,  there  is  much  better  agreement  with 
the  present  numerical  theory  than  with  Leehey's  asyii^>totic 
theory.  The  foil  and  cavity  combination  is  no  longer  of 
large  aspect  ratio  at  small  O'/o*^  » so  the  asymptotic 

theory  is  not  valid. 

The  mcanent  coefficient  is  taken  about  the  mid-chord, 
consistent  with  the  right  hand  rule.  Figures  4.6  and  4.7 
show  the  moment  coefficient  for  the  aspect  ratio  5 and  3 
hydrofoils.  Since  Leehey's  matched  asymptotic  expansion 
theory  neglects  the  lifting  surface  effect,  it  is  expected 
that  the  present  numerical  lifting  surface  theory  would 
show  better  agreement  with  experiments. 

The  nondimensional  cavity  length  versus  the  ratio  of 
cavitation  number  to  angle  of  attack  is  shown  in  Figures  4.8 
and  4.9.  The  cavity  length  is  measured  from  the  leading 
edge  at  the  spanwise  location  of  the  centroid  of  the  foil 
area.  Generally  speaking,  the  agreement  between  the  theory 
and  experiments  is  very  good. 

In  Maixner's  experiment,  the  cavity  pressure  was 
measured  with  a foil  surface  pressure  tap.  A noticeable 
"hook"  was  found  in  the  lift  and  moment  data  at  higher  emgles 
of  attack  (Figures  4.4  & 4.6  ).  Reun  effects  on  the  cavity 
pressure  measurement,  due  to  the  dyn2unic  pressure,  were 
further  investigated  on  the  large  foil  used  in  Maixner's 
test.  The  cavity  pressure  readings  are  taken  both  inside  the 
cavity,  using  a total-head  tube,  and  on  the  foil  surface.  The 
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L- shaped  total  head  tube  protuded  downwards  into  the  cavity 
from  the  upper  tunnel  wall  so  that  it  was  parallel  to  the 
foil  surface  and  pointed  towards  the  leading  edge,  away 
from  the  impinging  re-entrant  jet.  Figure  4.10  shows  that 
the  readings  from  the  toil  surface  pressure  tap  are  consis- 
tently higher  than  the  measured  cavity  pressure,  especially 
at  higher  cuigles  of  attack  and  shorter  cavity  lengths.  If 
the  correct  cavity  pressure  readings  were  taken  in  the 
experiments  of  Maixner,  the  discrepancy  from  theory  at  higher 
angles  of  attack  and  shorter  cavities  would  probeibly  have 
been  smaller. 

The  numerical  results  for  an  unsteady  two-dimensional 
hydrofoil  with  a finite  cavity  are  shown  in  Figures  4.11  to 
4.14.  Lift  and  moment  coefficients  in  heave  and  pitch  about 
the  leading  edge  are  compared  with  the  results  predicted  by 
Woods  (1957).  The  cavitation  number  was  zero  in  Wood's 
calculation.  At  zero  cavitation  number,  the  cavity  length 
change  was  neglected,  as  was  the  cavity  volume.  However, 
present  two-dimensional  solutions  are  the  inner  solution  of 
the  three-dimensional  case.  The  following  definitions  are 
used  for  the  unsteady  force  coefficients. 

Lift  due  to  Keo-ve  oLoa'I'  LB 
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The  phase  angles  are  defined  as  the  arc  tangent  of  the  imaginary 
part  to  the  real  part  in  the  last  equations. 

The  effect  of  finite  cavities  tends  to  reduce  the 
cunplitude  of  the  oscillating  force  and  moment  coefficients. 

The  diversity  of  phase  angles  between  finite  cavity  and 
infinite  cavity  increases  with  decreasing  reduced  frequency 
{‘dC/U»).  One  can  expect  this  happens  due  to  the  variation 
of  cavity  length,  it  increases  as  reduced  frequency  decreases. 

At  large  cavity  lengths  (low  cavitation  numbers) , 
the  oscillating  lift  force  was  away  from  the  region  of 
resoncuit  behavior.  Experimental  data  cuid  theoretical  values 
are  compared  with  present  numerical  results  in  Figures  4.15- 
4.18.  The  theoretical  values  for  cavitating  flow  over 
rectangular  foils  of  aspect  ratio  1 and  6 were  tedcen  from 
Widnall  (1966) , by  forming  a linear  comvination  of  heave  and 
pitch  about  the  leading  edge,  resulting  in  pure  pitch  about 


the  midchord.  The  maximum  chord  of  the  semi-eliptic  foil  was 
taken  as  the  chord  of  the  rectangular  foil.  The  moment 
coefficient  was  based  upon  planform  area  and  mean  chord. 

Figures  4.15  and  4.16  show  the  lift  coefficient  versus  reduced 
frequency  for  the  aspect  ratio  5 and  3 hydrofoils, 

respectively.  The  same  conditions  as  in  experiments  were 
used  in  the  calculations.  The  agreement  between  numerical 
values  and  experimental  results  is  good  for  lift  amplitude. 
However,  the  phase  eingle  predicted  by  theory  is  consistently 
higher  than  experiment.  This  may  be  in  part  due  to  the 
variation  of  cavity  pressure,  which  is  sensitive  to  the  phase 
angle.  Although  the  cavity  pressure  variations  were  small  at 
long  cavities,  the  cavity  length  is  very  sensitive  to  the 
cavitation  number  (Figures  4.8  and  4.9).  However,  the  explan- 
ation of  phase  angle  discrepancy  is  not  fully  satisfactory  at 
this  moment.  The  difference  of  unsteady  lift  measurements 
between  these  two  hydrofoils  is  negligible.  The  finite  span 
effect  is  not  noticable  for  the  frequency  rauvge  of  the  tests. 
Figures  4.17  and  4.18  show  the  oscillatory  moment  coefficient, 
for  the  aspect  ratio  5 and  3 hydrofoils.  At  high  reduced 
frequencies,  the  present  theory  predicts  larger  moment  coeffic- 
ients that  found  by  experiment.  The  difference  of  phase 
angles  between  theory  and  experiment  is  large.  One  possible 
reason  is  that  the  oscillatory  moment  of  supervacitating  flow 
is  a small  quantity,  which  has  the  same  order  of  magnitude  as 
the  moment  produced  by  the  inertia  force  and  spring  force  in 
the  mechanical  system.  The  torsional  system  of  the  dynauno- 
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meter  used  in  measuring  the  moments  was  not  calibrated  at 
the  conditions  that  the  unsteady  measurements  were  made.  The 
con^arison  between  the  present  numerical  method  for  a rectang- 
ular hydrofoil  and  Widnall's  results  is  shown  in  Appendix  E. 

Present  numerical  calculation  also  gives  the  variation 
of  cavity  length  and  cavity  volume.  Figures  4.19  and  4.20 
ahow  these  two  variations  decrease  with  reduced  frequency 
increases.  The  eunplitude  of  the  oscillatory  cavity  volume  is 
ckbout  40%  of  the  mean  cavity  volume  when  Aa/a  = .167.  If  the 
variations  of  cavity  pressure  were  considered,  the  phase  lag 
should  enlarge  due  to  the  heat  transfer  process. 
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Chapter  V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Measurements  were  made  of  forces  and  moments  for 
aspect  ratio  3 and  5 supercavitating  hydrofoils  of  ellip- 
tical planform  oscillating  in  pitch.  For  short  cavities, 
the  natural  frequency  of  the  cavity  plays  an  important  role 
in  the  oscillating  forces.  A pressure  transducer  was  used 
to  measure  the  variation  of  cavity  pressure.  It  showed 
that  the  cavity  pressure  is  not  constant  during  the  pitch 
oscillation.  The  amplitude  of  cavity  pressure  oscillation 
depends  on  the  cavity  length  and  oscillation  frequency. 

The  thermodynamic  analysis  of  the  cavity  was  performed 
by  assuming  spherical  symmetry  and  uniform  cavity  interior 
properties.  Eigenvalues  were  calculated  subject  to  the 
given  boundary  conditions.  The  natural  frequency  decreases 
when  the  cavity  volume  increases,  which  confirmed  the 
experimental  results.  The  stiffness  of  the  cavity  results 
from  the  air  content  inside  the  cavity.  The  calculated 
ratio  of  imaginary  part  of  the  natural  frequency  to  the 
undamped  natural  frequency  indicated  that  more  air  content 
decreases  the  damping  ratio.  The  temperature  gradient 
inside  the  cavity  was  neglected  in  the  analysis.  The 
model  used  only  describes  general  features  of  the  process. 
The  entire  unsteady  analysis  must  be  redone  to  account  for 
cavity  pressure  variations. 

The  discrete  vortex  and  source  method  is  developed  for 
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supercavitating  hydrofoils.  The  cavity  length  was  iterated 
to  get  a uniform  cavitation  number  over  the  cavity  planform. 
The  numerical  results  for  an  unsteady  two-dimensional 
hydrofoil  with  a finite  cavity  were  compared  with  the 
theoretical  results  at  zero  cavitation  number.  The 
appreciable  difference  of  phase  angles  between  non-zero 
and  zero  cavitation  numbers  was  found.  The  lift  and 
moment  calculations  for  supercavitating  hydrofoils  of 
elliptic  planform  in  steady  flow  was  performed  and  compared 
well  with  previous  experiments.  A more  accurate  prediction 
of  lift  and  moment  coefficients  was  obtained  by  the  present 
numerical  method  than  with  existing  asymptotic  theories. 

For  long  cavities  the  numerical  results  of  three-dimensional 
unsteady  flow  were  compared  with  experiment.  Amplitude 
of  oscillatory  forces  compared  well  between  theory  and 

I experiment.  However,  noticeable  difference  of  phase 

angles  was  found.  Since  the  cavity  pressure  was  considered 
as  a constant  during  unsteady  motion,  one  should  find  the 
importance  of  cavity  pressiire  variation. 
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Table 


2-D  Steady 


3-D  Steady 


2-D  Unsteady 


3-D  Unsteady 


4.1  Values  of  Vortex  and  Source  Element 


No.  of  vortex 
element 

No.  of  source 
element 

25 

40 

No.  of  element 
along  the  chord 
on  the  foil 

No.  of  element 
in  the  cavity 
behind  the  foil 

No.  of  element 
along  the 
semispan 

12 

8 to  14 

6 

No.  of  vortex 
element 

No.  of  source 
element 

20 

35  to  40 

No  of  element 
along  the  chord 
on  the  foil 

No.  of  element 
in  the  cavity 
behind  the  foil 

No.  of  element 
along  the 
semispan 

10 

8 to  14 

4 

Figure  2.1  Tunnel  Test  Section  and  Driving  System 


Figure  2.4  Instrumentation  Arrangement  for  Hydrofoil  Test 
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FIGURE  2.6  a Non-I 

Davitating  Oscillatory  Moment  Coefficient 
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FIGURE  2.8  a Oscillatory  Lift  Coefficient, 
0 = 12“  Aas±2®  AR=3 
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FIGURE  2.10  Oscillatory  Moment  Coefficient, 
a = 12®  Aa  = ± 2®  /R  = 3 
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FIGURE  2.10  Oscillatory  Moment  Coefficient, 
a « 12*  Aa  = ± 2*  /R  = 3 
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FIGURE  4.1  FOIL  PLANFORM  AND  TYPICAL  VORTEX  SOURCE 
ELEMENT 
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FIGURE  4.4  C^_/a  vs.  cr/a, 
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APPENDIX  A 

DYNAMIC  CALIBRATION  OF  DYNAMOMETER 

Dynamic  calibrations  of  the  water  tunnel  force  gages 
were  carried  out  with  the  aspect  ratio  5 foil.  The  following 
equipment  was  used  to  perform  the  calibration: 

1)  Vibration  Exciter  Body  - Bruel  & Kjaer  Type  4801, 

2)  Vibration  Exciter  Head  - Bruel  & Kjaer  Mode  Study 

Head  Type  4814, 

3)  Power  Amplifier  - Bruel  & Kjaer  Type  2707, 

4)  Impedance  Head  - Wilcoxon  Z-602  #252, 

5)  Sweep  Oscillator  - Spectral  Dynamics  SD  104A-5, 

6)  Dynamic  Analyzer  - Spectral  Dynamics  SD  101-A, 

7)  Buffer  Attenuator  - Home  Made, 

8)  Ithaco  Preamplifier  - Model  255-A, 

9)  Accelerometer  - Bruel  & Kjaer  Type  4332  #159861, 

10)  Phase  Meter  - Spectral  Dynamics  SD  110, 

11)  Oscilloscope  - Te)ctronix  Type  561B, 

12)  RMS  Voltmeter  - Bruel  & Kjaer  Type  2416, 

13)  DC  Voltmeter  - HP  Type  405  CR, 

14)  Transducer  Digital  Indicator  - Lebow  Model  7521. 

The  instrumentation  arrangement  is  shown  schematic- 
ally in  Fig.  A.l 


3 


. u»«i 


The  dynamometer,  installed  in  the  water  tunnel,  was 
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sub  jectied  to  a sinusoidally  oscillating  force,  generated 
from  B & K Exciter  Head  Type  4814.  The  exciter  was  mounted 

! 

on  a heavy  table  and  was  also  fixed  to  the  side  flange  of  | 

the  tunnel.  The  stiffness  in  the  connections  between  the  j 

exciter  2md  the  flange  of  the  tunnel  was  sufficiently  high 
that  the  resonance  frequency  was  out  of  the  range  of  interest. 

A pushrod  was  used  to  transfer  the  force  from  exciter  head 
to  the  pressure  center  of  the  foil:  an  impedance  head  was 
attached  between  them.  Figure  A. 2 shows  a picture  of  the 
exciter  arrcingement.  The  exciter  was  driven  by  the 
Oscillator  and  B & K Power  Amplifier,  which  is  a high  power, 
low  distortion  vibration  exciter  driver,  giving  a highly 
linear  output  and  a selectable  output  impedance. 

The  force  signals  from  the  load  cells  and  the  impedance 
head  were  analyzed  with  the  Spectral  Dynamics  Analyzer,  in 
which  a 5 Hz  filter  was  used.  The  signal  being  analyzed 
was  compared  internally  to  a signal  of  the  same  frequency 
as  the  exciting  frequency  from  the  oscillator  and  filtered 
with  the  5 Hz  Band  Pass  filter  around  that  frequency.  The 
measurement  covered  the  range  from  5 Hz  up  to  30  Hz.  Due 
to  problems  of  stabilization  of  the  shaker  output,  the  tests 
were  carried  out  for  discrete  frequencies. 

The  phase  cingle  between  the  force  outputs  of  the  load 
cell  amd  the  impedance  head  was  measured.  Both  of  these 
signals  come  through  the  Spectral  Dyneunics  Analyzer  and  were 
compared  in  the  Phase  Meter.  The  three  load  cells  on  the 


steel  rods  that  sustain  the  floating  section  of  the 
dynamometer  were  replaced  by  dummy  load  cells  to  increase 
the  stiffness. 


The  dynamic  calibration  of  the  torque  transducer  was 
the  most  difficult  calibration  made.  It  was  hard  to  calibrate 
the  moment  gage  under  the  same  conditions  that  the  unsteady 
measurements  were  made.  During  the  calibration,  the  top 
part  of  the  torque  transducer  was  rigidly  cleimped  to  the  base 
of  the  dynamometer  and  the  top  window  of  the  tunnel.  The 
moment  input  was  obtained  from  the  force  on  the  impedance 
head  multiplied  by  the  moment  arm.  The  output  of  the  torque 
transducer  was  obtained  from  the  Lebow  Indicator.  The 
calibration  showed  that  the  moment  from  torque  transducer 
had  a linear  relation  to  the  applied  moment. 

The  dynamic  system  of  the  dynamometer  in  lift  and  drag 
directions  can  be  described  by  a model  of  three  degrees  of 
freedom.  The  applicable  equations  of  motion  are 

VI, a,  + + CA,  -Xj)  t Xj')  -D , 

, . _ j«#t 

THiXi  + "ft*  , 
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where 

= mass  of  the  base  of  the  dynamometer 
and  top  window  of  the  tunnel » 
m^^  = mass  of  the  floating  part  of  the.  dyncunometer 
and  23%  of  the  shaft  mass, 
m2  = mass  of  the  foil  and  77%  of  the  shaft  mass, 
k2,C3  = spring  constant  and  damping  coefficient 
in  the  tvmnel's  window,  respectively, 
kj^,Cj^  = spring  constant  and  damping  coefficient  of 
the  load  cell  plus  the  linkage  between  the 
load  cell  and  the  dynaimoneter , respectively, 
k2,C2  = spring  constant  and  damping  coefficient  of 

the  shaft  and  foil  combination,  respectively. 

The  shaft  mass  was  divided  as  described  in  Den  Hartog 
(1956)  and  based  on  Rayleigh's  method. 

A contilever  beaun  El,  of  length  jl  and  mass  m carries  a 
concentrated  mass  M at  its  end.  The  deformation  of  the  beam 
is  approximated  by  a quarter  cosine 


The  total  kinetic  energy  cam  be  written  as 


K.6.  = j . 
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The  potential  energy  will  be 

K.E.  = (H  + 0.22.7/)). 

Equating  the  kinetic  energy  to  the  potential  energy  and 
solving  for  oo^  , we  find 

.2  3.o3£I 

U)  “ —z 

Thus  23  percent  of  the  mass  of  the  beam  is  effectively 
added  to  the  end  mass. 

Equations  (A-1)  can  be  solved  by  assuming  harmonic 
motion  for  the  displacements.  The  ratio  of  Lebow  force  to 
exciter  force  can  be  obtained  after  the  spring  constants 
cuid  masses  are  known.  The  value  of  these  constants  are; 

Mass  of  the  base  of  dynamometer  and  top  window  of 
the  tunnel 
m^  = 10  slug. 

Mass  of  the  floating  section  of  dynamometer  = 

m^  = 1.2  slug. 

Mass  of  the  shaft 

m =1.37  lb  =0.043  slug, 
s m 

Mass  of  the  floating  part  of  the  dyneunometer  and 

a fraction  of  the  shaft  mass 

m.  = m-  + 0.77  m = 

1 r s 

1.23  slug f 


1.2  + 0.043  X 0.76 
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Mass  of  the  hydrofoil 

= 0.157  slug. 

Mass  of  the  foil  and  a fraction  of  the  shaft  mass 
m2  = mj^  + 0.23  X m^  = 0.168  slug. 

The  load  cells  are  coupled  between  the  floating  part 
and  the  base  of  the  dynamometer  by  screw  rods  eind  bolts. 

The  spring  constant  of  the  load  cells  was  obtained  from  the 
Lebow  catalog.  The  spring  constant  of  the  rods  are  calcu- 
lated from  Hooke's  law,  neglecting  the  Poisson  effect.  The 
equaivalent  spring  constant  of  the  lift  load  cell  system 
(#125)  is 

= 66440  Ib/in. 

The  spring  constant  of  the  base  dynamometer  and  top 
window  combination,  was  obtained  from  the  experimental 
force-displacement  curve.  This  work  was  done  by  W.  T.  Tan 
as  part  of  UROP  project.  This  constant  is 

k^  = 299,000  Ib/in. 

Spring  constant  of  the  shetft  and  foil  is  (Riva  (1974)) 
kj  = 25,000  Ib/in. 

Damping  is  small  in  this  mechanical  system.  1%  of 
the  individual  critical  damping  of  each  spring-mass  system 


is 
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= 1.65  Ib^-sec/in, 

C2  = 0.375  Ib^  sec/in, 

= 10  Ib^  sec/in . 

From  the  calibration,  it  is  evident  that  the 
natural  frequency  of  the  dynamometer  is  above  the  test 
range.  The  results  show  that  it  is  reasonable  to  represent 
the  real  system  by  the  model  with  three  degrees  of  freedom. 

It  also  shows  that  the  damping  is  small  in  the  system  and 
the  dyncimometer  can  be  modeled  as  a mass-spring  only.  The 
lowest  natural  frequency  of  the  dynamometer  with  foil  a 
attached  was  about  80  Hz.  The  maximum  frequency  of  the 
pitching  oscillation  was  about  20  Hz.  The  mass  of  the  foil 
system  and  of  portions  of  the  dynamometer  produced  a force 
on  the  load  cells  when  the  apparatus  was  oscillated.  The 
magnitude  of  the  force  was  determined  by  oscillating  the 
foil  in  air  and  measuring  the  resulting  force.  This  force 
was  sTobtracted  from  the  response  of  the  supercavitating 
foil,  taking  into  consideration  the  phase  relationship 
between  the  two  forces. 

The  torsional  system  can  be  represented  by  a model  with 


two  degrees  of  freedom.  Equations  of  motion  are 
D",  ^ ^ t ^^5  - Vi)  + Av  V',  = M, , 


in  K M, 


fA-2) 
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where  and  J2  are  the  moment  of  inertia  of  the  foil  and 
shaft  combination  and  crankarm  respectively;  and  M2  are 
the  moment  applied  at  the  foil  and  crankarm,  respectively; 
c^  amd  k^  are  the  damping  coefficient  and  torsional 
spring  consteint  due  to  the  bearing  and  rubber  sealing.  The 
moment  meausred  from  the  torque  trandsucer  will  be 


Substituting  this  into  the  last  set  of  equations.  Using 

*9  •• 

the  assumption  of  simple  harmonic  motion  and  Tf,  - Vi  , these 
equations  become 


C^'3) 


where  the  phase  angles  » are  relative  to  the  acceler- 

ation. These  equations  can  be  simplified  as 

'Jt ^ V;<. e”*** - +1,. H. e ^ * 


If  the  constamts  Jj^,  and  k^  are  known,  then  the  moment 
applied  at  the  foil  can  be  calculated  from  the  equation 
(A-4)  after  the  torque  transducer  output  is  measured. 


-106- 


When  the  foil  was  oscillating  in  the  air,  the  moment 
applied  to  the  foil,  was  zero.  The  equations  of  motion 
become 


dt" 


= 0, 

4 Mt  = M,. 


r/>-5-; 


The  moment  M2  was  measured  by  the  force  at  the  load  cell 
in  the  drag  direction  multiplied  by  the  length  of  the 
crankarm.  was  measured  by  the  torque  transducer.  The 

phase  relationship  was  measured  with  respect  to  the  acceler- 
ation. These  equations  can  determine  the  characteristic 
constants  of  the  system.  The  results  of  the  calibration 
show  the  damping  coefficient  to  be  dependent  on  the 
frequency.  For  the  aspect  ratio  3 hydrofoil,  these  constants 
are 


O',  0.00  24  -ff-lk'Src', 

= 0.0033  -ff-lW-Sec*, 

=4.^3  ^+-lb. 


(A- 6) 


From  the  results  of  the  calibration,  it  is  reasonable  to 
represent  the  torsional  system  by  a model  with  two  degrees 


of  freedom. 


Figure  A. 1 Instrumation  Arrangement  for  Calibration 


[ 

I 

I 


I 

j 

I 

t 

f 

k 

I 


f 


L 


-108- 

APPENDIX  B 

FORMULATION  OF  PROBLEM 


In  the  following,  the  isentropic  inviscid  flow  of 
a incompressible,  initially  irrotational  fluid  is  considered 
for  linearized  supercavitating  hydrofoils.  Under  this 
hypothesis,  the  flow  can  be  described  by  the  velocity  potent- 
ial $ . The  linear  equation  of  the  velocity  potential 
exterior  to  the  foil  plus  the  cavity  plus  the  vortex  sheet 
is 

(B-1) 

where  V*  is  the  Laplacian  operator.  It  is  convenient  to 
introduce  the  perturbation  potential  y , such  that 

+ (3-2) 

where  u.  is  the  undisturbed  flow  velocity  in  the  direction 
of  the  positive  x-axis.  Combining  equations  (B-1)  and  (B-2) 
yeilds  the  equation  for  the  perturbation  potential 

vV  = o.  (3-3) 

The  linearized  Bernoulli  equation  for  unsteady  flow 


IS 


f-r.  --?(■ 


at 


(3-4) 


which  yields,  for  the  pressure  coefficient 

-t  ~ f -L -I- ) (5-5) 

The  perturbation  pressure,  P , also  satisfies  the  Laplace 


equation 


V*P  = 0, 
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r 


exterior  to  the  foil  plus  the  cavity  plus  the  vortex  sheet. 

Let  the  source  distribution,  » represent  the  slope 

difference  between  the  upper  and  lower  surface  at  point 
and  the  vortex  distribution , represent  the 

difference  in  the  x-component  of  perturbation  velocity  between 
the  upper  and  lower  surface  at  point  As  shown 

in  Fig.  B-1,  this  leads  to  the  results, 

o S (S-7) 

^ 

(b-8) 

Oa''’''0.tJ  ' U«  ' » 

where  denotes  the  jump  in  "a"  across  the  surface  ^*0. 

From  the  Bernoulli  equation  (B-5)  and  the  distributions 
of  sources  and  vortices,  one  can  see  that 


= (5-9) 

^ ^ U,  at  ax 

X '• 


9 

. J 

f* 


4 


i 


\ i 

L 
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For  the  application  of  Green's  theorem  to  the  super- 
cavitating  hydrofoil  of  finite  span,  the  pressure  at  the 
point  J)  is  (e.g.  Newman  (1977)) 

where  r is  the  distance  between  and  dummy  point 

and  ^ encompasses  cavity,  foil  and  wake. 

The  integral  can  be  linearized  for  application  to 
thin  wings  at  small  angles  of  attack  and  rewritten  as  open- 


surface  integral 


= ( <p>i  f -<|£>  , 


(6-/2) 


where  5^  is  the  projection  of  the  foil  surface  on  to  the 
X2,  pland  and  5^  is  the  projection  of  cavity  surface  on  to 
the  Scune  plane. 


The  partial  derivative  -ip-  , is  equal  to  zero  everywhere 
on  the  plane  J=o  except  at  the  point  . If  we  assume  the 

pressure  distribution  is  continuous  over  the  surface,  then 
the  limiting  case  for  r=o  can  be  solved  as  in  the  following. 


a-*  x-f 


f3-/3) 


} 

i 

{ 
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If  the  point  is  on  the  surface  ^ ^ O # by  introducing 
temporary  integration  variedjles  "5' = fx- "5  ),  ■5'=  ( J-5)  , we  find 


~i  ■* 


The  limiting  process  involved  here  shows  that  we  must  let 
^ and  6 approach  zero  in  such  a way  that  the  ratio 
becomes  indifinitely  large  which  leads  to 


(5-/3) 


Substituting  equations  (B-9) , (B-10)  and  (B-15)  into  Eq. 
{B-12) , the  pressure  on  the  surface  becomes 

Sis.  * 

•« 

The  second  tem  in  the  first  integral  can  be  integrated 
by  parts 


I] 


- J , (6-7) 


^iU£'  .lL!*Ui4JU|ilifa.l-9ia0V 


-112- 


where  we  use  the  condition  that  the  source  distribution 
is  zero  outside  the  cavity. 

When  the  boundary  condition  of  constant  pressure 
in  the  cavity  is  applied  on  the  cavity  surface,  the  follow- 
ing result  is  obtained 

(e>-IB) 


i fU." 


By  combining  Eqs.  (B-16) , (B-17)  and  (B-18) , the  integral 
equation  becomes 


+ 2^ 


J 


j I p o-at  rJJ5. 


(6-/9) 


This  integral  equation  can  also  be  derived  by  using 
the  vortex  potential,  source  potential  and  Bernoulli  equation. 


Applying  the  Green's  theorem  to  the  Laplace  equation 
, the  velocity  potential  at  an  arbitary  field 
point  is 

7)^5- 
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Again,  the  surface  integral  encompasses  the  cavity,  foil 
and  wake.  Since  the  operation  of  differentiation  with 
respect  to  k can  be  interchcinged  with  V*  in  Eq.  (B-3) , 
we  may  write 


-f^-p)dS.  (8-2/) 


In  regard  to  the  second  term  in  Eq.  (B-21) , the 


condition  of  irrotationality  reveals  that 

, 


(fi-22) 


on  top  and 


au  _ 

; 


C5-a3) 


on  the  bottom  surface.  Eq.  (B-20)  can  be  rewritten  as 
open  surface  integral  for  linearized  problem, 

jj  (<“>^7-<ln>'rl 

- ± ff  ( r ^ -i-  -il  d5  (B-2^) 

"4n  11  r*«37  r ^9%'  r}‘‘  ■ 

Si+5, 

From  the  relation  U =0(P/‘3X  and  V — 


(3-2S) 
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where  the  assumption  that  ^ = 0 has  been 

made  without  loss  of  generality.  Inserting  Eg.  (B-25)  into 
Eq.  (B-24)  and  interchanging  orders  of  differentiation  emd 
integration,  we  get 


4tt 


J 


(b-26) 


The  integration  of  first  term  is  zero  everywhere  on  the  plane 
^ = 0 except  at  the  point  f - O . The  limiting  value  of 
this  has  been  shown  in  Eq.  (B-15) . The  result  of  the 
second  term  in  the  integral  is  the  Mangier' s principal 
value  which  is  equal  to  (Ashley  & Landall  (1965)) 
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By  combining  Eqs.  (B-7), (B-15)  and  (B-27) , the  integral 
equation  which  satisfies  the  tangential  velocity  on  the 
boundary  becomes 
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For  the  calculation  of  oscillatory  forces  and  moments 
due  to  small  amplitude  oscillation  about  the  steady  state, 
the  integral  equations  (B-19) , (B-28)  are  solved  for  subject 
to  the  closure  condition 

ds  = 0.  (B-29) 

A Kutta  condition  is  also  required  to  make  the 
problem  well  set. 

The  unknowns  are  the  vortex  strength,  source  strength 
and  cavity  surface. 

As  mentioned  before,  the  cavitating  lifting  surface 
integral  equations  (B-19)  and  (B-28)  can  also  be  derived 
by  using  vortex  and  source  potentials.  A sheet  of  source 
spread  continuously  over  the  surface  Sc  , having  strength 

per  unit  area  in  the  neighborhood  of  the  point  , 

possesses  the  perturbation  potential 


(S-30) 


The  vertical  velocity  at  point  ,-o,  3.)  induced  by 
the  source  sheet  is 
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The  pressure  at  point  (^>  + 0.^)  induced  by  the  source 
sheet  can  be  calculated  from  Bernoulli  equation 


n -r. 

i eu; 


= - -2_  / J_  2!&. 


+ 


\ 


= -2 


[± 

(Ll 


(-4S 


4Tr  dt 


JJ 


p'i.9.±n'3d5  1 


( 

1 :§_ 

4it9x 

J. 

^ 

Uo  9-t  /V  <'<5-:?;"]'^ 


(&-32) 


-117- 


Instead  of  solving  the  velocity  potential  directly, 
we  may  distribute  the  vortices  over  the  foil  and  the 
weUce , Sw,  to  form  the  integral  equation.  The  vertical 
velocities  induced  by  a vortex  sheet  are  determined  by 
using  the  Biot-Savart  law,  (Bisplinghoff , et  al  (1957)) 
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where  JT  and  £ are  bound  vortex  strength  and  trailing 
vortex  strength,  respectively.  The  Kirchhoff's  law  amd 
Kelvin's  theorem  must  be  satisfied  for  the  vortex  distri- 
bution. The  pressure  due  to  the  vortex  distribution  is 

p =.  _ 2.  r-L  ^ ) (3-34) 

where  d^v/dX  is  the  x - component  velocity  induced  by  the 
vortex  distribution 


(3-35) 


(5-3?) 


This  is  equivalent  to  Eq.  (B-28) . The  same  expression  will 
be  obtained  if  we  apply  the  Kirchhoff's  law  and  Kelvin's 
Theorem  over  the  surface  and  evaluate  the  integration.  The 


of  uniting  Eqs.  (B-32)  and  (B-38) 


which  is  the  Scune  as 


is  also  time  dependent 


Note  that  the  cavity  region 
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APPENDIX  C 

VELOCITY  INDUCED  BY  A DISCRETE  VORTEX  SEGMENT 

Consider  the  vortex  segment  with 

strength  JT  per  unit  length.  Defining  as  the 

coordinates  of  a general  point  on  the  vortex,  and  Cx,jr)  as  the 
coordinates  of  a control  point.  The  law  of  Biot-Savart  is 
applied  to  obtain  the  induced  velocity  at  point 

Along  the  vortex,  we  have 

-t  s = Cor>s+an+, 

and  the  geometrical  relationships 


a * 1 + 

e = + , 

= (x-x,+id,-f5)d5 -C3-f)-tdS 
= d-di  , 

[u-'? ( j - ? = Id  I /^>4  sec<^  , 

AS  = \d\/a.  . 
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Then  the  equation  (C-1)  becomes 

9* 

= — ( , <'C-3) 

where  o<  and  ^ are  the  angles  shown  in  Fig.  (C  .1) 

Equation  (C-3)  is  not  suitable  for  numerical  ccsnputa- 
tion  if  d becomes  small  (Kerwin  & Oppenheim  (1974)). 

There  are  three  special  cases  which  require  numerical  approx- 
imations for  small  cl  , 
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ccycX  = - CO30  ^ ~l  + 


(3)  for  ^ > S'! 

oL  ^ 1 2d/(2.ai,-tb^  1 , 0 ^TT- f * t2({/(i«  Ji+k')!  , 

ccsoi  ^ 1 - adVfaaa.+  bf  , 

COS^  = -£OS0  = -I  + 2(^yC2fl^a  + b)^, 


4^  V 


= 2.a'^\ z-  T-^ il  ^C-6) 


The  induced  velocity  to  the  trailing  vortex  iTii.O.i,") 
“C^J/OjJO  with  strength  Sa  is 


41^^  = - 


^ 7<-^  ! 

3"  a. 

t-<. 

The  total  velocity  induced  by  a vortex  distribution 
is  the  sum  of  all  the  spanwise  and  chordwise  vortices. 
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FIGURE  C.l  Coordinate  Systems  for  Discrete 
Vortex  Segment  and  Control  Point 
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APPENDIX  D 

PRESSURE  INDUCED  BY  A DISCRETE  SOURCE  SEGMENT 

For  unsteady  supercavity  flow,  the  pressure  Induced 
by  source  component  consists  of  two  terms:  the  space 
derivative  3</k/ax  and  time  derivative  aVi/at- . Rewrite 
the  Bernoulli  equation 


p Ts  ~ Ta  _ ^ f f \ 


The  velocity  in  the  x - direction  at  C’f.j)  due  to  the 
line  source,  ^ (t)  , along  CX, , ) - f if,,  j,')  is  (Fig.D.l) 
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where 
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and  (t,b.  d 


are  defined  the  same  as  Eq.  (C-2) 
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Again,  there  are  three  special  cases  which  require 
numerical  approximation  for  small  d . Equation  (D-1)  can 
be  represented  by  the  geometric  emgles  (Fig.  D.2) 
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(3)  for  }f>2tz 

COSC^~  I - 2rfyC2Gj,-fbf  , S:nol  =i\2d/(:>.al,^hi\  ^ 
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The  induced  pressure  due  to  time  variation  of  the 
source  segment,  , is 
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By  using  the  definition  of  Eq.  (C-2) , we  obtain 
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Once  more,  three  special  formulas  are  considered  for 
small  d I 
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(3)  for  ^ 

<I>I=  -(f 
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Then  the  induced  pressure  due  to  source  distribution 
can  be  calculated  by  using  these  equations  and  Bernoulli's 
equation. 
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APPENDIX  E 

NUMERICAL  RESULTS  OF  RECTANGULAR  HYDROFOIL 

The  present  nxunerical  method  was  used  to  calculate  the 
unsteady  forces  on  an  aspect  ratio  6,  rectangular  hydrofoil. 

A comparison  was  made  on  heaving  motion  at  a /a  = 2.0  and 
wCq/U_^  = 3.1416.  Eight  elements  were  used  along  the  chord 
of  the  foil,  six  elements  were  used  in  the  cavity  behind  the 
foil  and  two  elements  were  used  along  the  semispan  in  the 
calculation.  The  lift  eimplitude  eind  phase  angle  were  calcu- 
lated after  one  cycle  of  oscillation  and  were  found  to  be  2.71 
euid  108.6°,  respectively.  The  lift  amplitude  emd  phase  angle 
of  Widnall's  results  at  wC  /U  = 3.1416  are  3.4  and  115°. 

O 00 


\ 


Ashley,  H.  & Landahl,  M.  1965  Aerodynamics  of  wings  and  bodies 
Addison-Wesley . 


H.  & Halfman,  R.  L.  1957  Aero 


Den  Hartog,  J.  P.  1956  Mechanical  vibration.  McGraw-Hill 


Kerwin,  J.  E.  & Oppenheim,  B.  W.  1974  A lifting-surface 

program  for  trapezoidal  control  surfaces  with  flaps 
Rep.  No.  74-15,  MIT 


Newman,  J.  N.  1977  Marine  hydrodynamics.  MIT  Press 


